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1.0  Project  objectives: 

An  important  part  of  the  design  process  for  turbomachinery  components  in  aeroengine 
applications  is  ensuring  that  the  component  is  resistant  to  failure  from  foreign  object  damage 
(FOD).  This  can  be  accomplished  in  a  cost-effective  manner  by  using  phyacs-based  models  in 
conjunction  with  experimental  testa  Thisconcept  isknown  as  Integrated  Computational  Material 
Sd  ence  and  Engi  neeri  ng  ( I  CM  SE) .  The  resul  ts  f  rom  anal  yti  c  model  s  are  used  to  gui  de  the  desi  gn 
of  a  test  matrix  and  thus  reduce  the  cost  of  the  experimental  testa  The  experimental  tests  are 
used  to  val  i  date  the  model  s  or  conversel  y  demonstrate  where  and  how  the  phya  cal  theory  needs 
to  be  redeveloped. 

Experimental  testing  can  be  expensive  and  time  consuming.  Conversely,  the  development  of 
model  sand  running  of  a  mulations  can  be  conducted  in  amatter  of  days  at  asmall  fraction  of  the 
cost  of  the  experiment  a  The  development  of  a  sound  design  requires  agreement  between  the 
theoretical  models  and  the  experimental  data  A  design  process  that  utilizes  experimental, 
theoretical  and  computational  tools  in  this  manner  is  the  only  practical  way  CMC  material 
systems  can  be  de/el  oped  and  used  in  high  performance  applications  in  a  cost-effective,  safe  and 
reliable  manner. 

The  objective  of  phase  I  was  to  develop  and  demonstrate  an  initial  methodology  for  analyzing 
the  effects  of  foragn  object  damage  (FOD)  on  aeroengine  components  utilizing  CMC  materials. 
The  methodology  integrates  the  effects  of  material  behavior  into  the  structural  design  process. 
The  long  term  objective  is  to  incorporate  this  analysis  capability  into  a  practical  engineering  tool 
for  ICSME  to  develop,  design  and  manufacture  the  components  for  tomorrow’s  high 
performance  vehicles.  This  methodology  is  an  enabling  technology  to  retan  a  national  technical 
capabi  I  i  ty  i  n  aerospace  engi  neeri  ng. 

2.0  Project  approach: 

The  primary  reasons  why  ceramics  have  not  been  widely  adopted  in  aeroengine  applications  are: 
(a)  limitations  in  tha'r  reliability  within  the  harsh  gas  turbine  environment,  (b)  interfacing 
ceramic  components  with  adjacent  metallic  components  leads  to  a  mismatch  of  material 
properties,  and  (c)  the  current  high  cost  of  manufacturing  components  New  ICMSE  tools  are 
required  to  overcome  these  obstacles.  A  key  element  of  this  particular  ICMSE  activity  is  the 
anal  ysi  s  methodol  ogy  used  to  model  the  FOD  and  i  ts  i  mpact  on  component  reliabtl  ity.  Thus  the 
phase  I  effort  is  focused  on  demonstrating  an  initial  capability  to  model  FOD  and  the  subsequent 
crack  and  fatigue  behaviors  in  CMC  components  The  first  step  in  developing  this  methodology 
was  to  re/iew  the  literature  on  FOD  to  identify  the  key  physical  phenomenon  present  in  FOD 
scenarios  The  objective  was  to  ensure  that  the  mechanisms,  which  need  to  be  engineered  in 
order  to  ensure  a  FOD  resistant  material,  are  properly  modeled.  The  approach  that  was  used  in 
phase  I  was  to  develop  physics-based  models  of  the  identified  constituent  material  and  fiber 
architectural  behavior  and  attempt  to  predict  the  damage  map  of  CMC  materials  The  results 
were  then  compared  to  available  experimental  data  to  evaluate  the  effectiveness  of  the  new 
model  sand  their  potential  for  further  development  into  an  ICMSE  tool. 

Figure  1  shows  the  framework  for  an  ICMSE  analysis  tool  that  can  be  used  for  evaluating  FOD. 
This  methodology  can  be  used  to  derive  a  potentially  optimal  CMC  materials  architecture. 
Which  will  allow  adesign  engineer  e/aluate  a  large  design  space  in  acost  effective  manner. 
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Figure  1:  Framework  of  analysis  codes  and  models  Starting  with  themodelsthat 
generate  the  I  ami  nate  properti  es  and  compl  eti  ng  wi  th  the  model  s  whi  ch  mi  mi  c  the 
physical  FOD  experiments 

The  design  process  must  define  the  internal  structure  or  material  architecture  for  the  CMC 
component  that  will  best  meet  the  structural  constraints  This  requires  that  the  physics-based 
predictive  tools  being  developed  can  accurately  predict  material  behavior  and  the  failure  modes 
due  to  FOD  in  composite  structures  However  to  date,  no  model  exists  to  describe  FOD 
phenomena  due  to  impact  dynamics  in  CMC  structures  This  can  be  attributed  to  the  complex 
nature  of  impact  dynamics  coupled  with  the  materials  architectural /constituent  complications 
Accurate  prediction  of  crack  behavior  includes  crack  initiation  and  propagation  for  both  single 
and  multiple  cracks  This  is  challenging  in  a  monolithic  material  alone,  and  many  approaches  to 
modeling  these  phenomena  require  prior  knowledge  of  crack  geometry  and  propagation.  This 
requires  additional  experimental  testing  to  acquire  this  knowledge,  where  a  truly  predictive 
capability  could  avoid  this  step.  In  a  composite  structure  the  situation  is  complicated  due  to  the 
potential  for  numerous  interfaces  between  fibers  fiber  coatings  and  lamina  Each  can  have 
different  material  properties  and  interactions  Key  failure  modes  that  must  be  resolved  in  a 
composite  structural  analysis  include  delamination,  fiber/matrix  debonding,  fiber  breakage,  and 
matrix  cracking. 

The  F’eridynamic  (F*D)  methodology  [Si Mings  2000]  was  evaluated  as  the  basic  methodology  for 
the  i  mpact  s  mul  ati  on  due  to  i  ts  uni  que  capabi  I  i  ti  es  to  meet  the  requi  rements  I  i  sted  above  T wo 
existing  computer  codes  were  investigated  to  provide  FID  with  laminate  material  properties: 
pcGINA  and  MAC/GMC.  The  integration  of  these  3  tools  (FID,  pcGINA,  MAC/GMC)  provides 
the  compl  ete  f  ramework  to  conduct  the  anal  ysi  s  as  descri  bed  above 
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3.0  Work  completed: 


3.1.0  Identification  of  key  physical  phenomenon: 

An  investigation  of  FOD  literature  has  been  conducted  to  identify  the  key  physical  phenomenon 
present  in  FOD  scenarios  The  objective  is  to  ensure  that  the  mechanisms  which  need  to  be 
engineered  in  order  to  ensure  a  FOD  resistant  material,  are  properly  modeled.  Thus  the 
developed  methodology  and  physic-based  models  need  the  capability  to  represent  these  phyg'cal 
behaviors 

One  caveat  is  in  order:  while  much  engineering  can  be  conducted  on  a  structural  level  to  meet 
the  life  constraints  of  the  material,  the  focus  of  this  work  is  to  develop  the  tods  which  will 
facilitate  the  tailoring  of  the  material  system  to  be  universally  applicable. 

To  engineer  FOD  resistant  materials  one  needs  models  which  are  able  to  simulate  the  proper 
damage  map  and  to  assess  the  residual  strength  of  the  material.  To  accomplish  this  task  the 
models  need  to  properly  represent  crack  formation  and  growth.  Proper  control  of  the  fiber/matrix 
interface  characteristic  isthe  key  to  understanding  crack  propagation  and  thus  obtaining  atough 
ceramic  composite 

3.1.1  Conical  damage  profile: 

In  [Chd  2008]  and  related  papers  FOD  was  experimentally  investigated.  Two  support 
conditions  were  considered:  first,  where  the  target  material  was  fully  support  and  secondly, 
where  the  target  material  was  only  partial  supported.  Furthermors  tests  were  conducted  at  both 
an  ambient  temperature  and  at  an  elevated  temperature  1316  C.  In  addition,  a  range  of  materials 
with  differing  hardness  were  cons'dered  for  the  impacting  object,  some  were  deformable  and 
others  were  brittle. 

Figure  2  shows  the  deta'ls  of  the  damage  resulting  from  the  impact  of  a  projectile  on  a  CMC 
materi  a I  sampl  a  There  i  s  damage  i  n  the  form  of  a  crater  at  the  i  impact  s' ta  Thi  s  damage  i  s  due 
to  the  breakage  of  matrix  and  fiber  tows  Furthermore,  cracks  with  cone  morphologies  form, 
which  emanate  at  approximate  120°  from  the  impact  crater.  Delamination  is  present  and 
spal  I  ati on  of  materi  al  from  the  backs' de  occurred  i  n  the  parti  al  supported  test  cases 

Figure  3  shows  an  idealization  of  the  damage  progress)  on  for  the  various  test  cases  considered  in 
[Choi  2008].  The  PD  results  demonstrate  these  patterns  very  well. 

[Choi  2008]  concluded  that:  impact  damage  increases  with  impact  velocity  and  that  impact 
damage  is  a  function  of  impactor  hardness  There  is  spallation  of  material  in  the  partial  support 
cases  due  to  an  additional  tensile  load  in  the  backside  region.  Lastly,  delamination  is  clearly 
evident  in  both  support  cases 
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Fi  gure  2:  Cross-sedi  onal  vi  av  of  the  i  mpact  si  te  of  a  S  C/S  C  target  sped  men  i  mpact ed 
at  400  m/s  by  1 .59-mm  steel  bal  I  proj  ecti  I  e  at  ambi  ent  temperature  Repn  nted  from 
[Choi  2008] 


[Ogi  et  al  2010]  conducted  similar  FOD  experiments,  but  they  also  considered  aCMC  with  a3D 
architecture  They  condudethat  tensile  stress  during  impact  was  transferred  to  thez-yarns  in  the 
3D-CMC,  which  constrain  del  ami  nation. 


Impact  site 


Backside  cracking 
&  delamination 


i 


Cone  cracks 


- A  / - 

- II — jZP - — — ~ 

/  V 

b=r— 

Low  velocity 


High  velocity 


(a)  Full  support  (b)  Partial  support 

Figure  3:  Schematics  of  damage  progress  from  low  (top)  to  high  (bottom)  impact  veloaties 
in  (a)  full  support  and  (b)  partial  support.  The  vertical  arrows  indicate  a  case  of  increasein 
impact  veloaty.  Reprinted  from  [Choi  2008] 


3.1.2  Mechanisms  of  fracture  toughness: 

In  metals,  plastic  deformations  occur  around  a  crack  tip,  which  allows  for  some  strain  energy 
dissipation.  Conversely,  monolithic  ceramics  have  brittle  failure  patterns  Thus  reinforcement 
f  i  bers  ar e  i  nd  uded  to  i  ntroduce  mechani  sms  f  or  the  di  ssi  pati  on  of  energy. 
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For  CMCs,  the  local  response  of  the  fiber  and  matrix  interface  is  critical  for  determining  the 
fracture  mechanics.  If  there  is  a  strong  bond  between  a  fiber  and  the  matrix,  then  a  crack 
propagating  through  the  matrix  has  a  bridge  or  access  point  into  the  fiber.  This  permits  the  fiber 
to  be  readily  cleaved  or  cut  without  atfecting  the  dynamics  of  the  crack  growth,  resulting  in  a 
lossof  load  carrying  capability.  Conversely,  if  this  interface  bond  is  weak,  then  the  propagating 
crack  might  cause  the  interface  to  break  and  the  fiber  and  matrix  to  debond.  This  allows  for 
several  energy-dissipation  phenomena  to  occur:  the  crack  may  be  deflected  by  the  fiber  and 
begin  to  run  perpendicular  to  its  original  path.  If  the  crack  does  not  change  direction,  then  it  will 
go  around  the  f  i ber  and  the  fiber  will  begi n  to  bri dge  the  crack  gap.  Dependi ng  on  the  nature  of 
the  crack  patterns  a  f  i  ber  may  debond  f  rom  the  matri  x  materi  al  at  several  poi  nts  al  ong  i  ts  I  ength. 
This  will  lead  to  a  situation  where  there  is  frictional  sliding  of  fiber  along  the  matrix  material. 
As  this  structured  cracking  of  the  matrix  material  occurs,  the  fibers  will  continue  to  carry 
increased  loads  This  will  lead  to  fiber  cuts  However,  in  this  case  the  amount  of  force  required  to 
break  the  fibers  will  be  greater  than  if  the  fiber  had  been  cleaved  due  a  strong  bond  with  the 
mdrix.  After  a  fiber  cut  has  occurred,  then  it  is  possible  for  the  fibers  to  pull  out  of  matrix.  As 
thi s  phenomena  occurs  one  will  also  seef  ricti onal  shear  stresses  devel op  at  the  i  nterf  ace 

3.1.3  Impact  resistance  of  EBC  coated  SC/SC  composites 

The  phase  one  analysis  did  not  explicitly  model  the  interaction  between  the  environment  barrier 
coating  (EBC)  and  the  CMC  substrate  material.  These  interactions  will  be  considered  in  a  phase 
2  analysis.  The  EBC  coating,  figure  4,  consists  of  three  sub  layers  first  a  bond  coat  layer  of 
silicon  is  deposited  on  top  of  the  substrate  followed  by  an  intermediate  mixed  layer  of  mullite  + 
bari  urn  stronti  urn  al  umi  num  si  I  i  cate  (BSA  S) ,  and  then  there  i  s  a  top  I  ayer  of  BSA  S. 


Figure4:  SEM  micrograph  of  typical  cross  section  of  plasma  sprayed  EBC  coating 
on  Ml  SC/SC  substrate  showing  microstructure,  composition  and  thickness  of 
EBC  sub-layers  Reprinted  from  [Bhatt  et  al  2008] 
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It  is  critical  to  raise  this  point  since  CMCs  require  reliable  EBC  coatings  to  operate  at 
temperatures  greeter  than  1 100  °C.  Without  an  EBC  coating  the  substrate  material  does  not  have 
adequate  life  [Bhatt  et  al  2008].  The  weakest  link  in  the  (EBC)  is  the  interface  between  the 
intermediate  coat  (a  mixture  of  mulliteand  BSAS)  and  the  silicon  bond  coating.  Any  cracking  or 
scratchi ng  off  of  the  EBC  provides  an  access  point  for  hostile  environmental  attack  to  the  CMC 
material.  It  is  worth  exploring  if  a  coating  can  be  used  to  increase  the  duration  of  impact  or  to 
reduce  impact  based  stresses. 

3.2.0  Material  selection 

The  material  chosen  for  analysis  is  a  melt-infiltrated  silicon  carbide  fiber/si  I  icon  carbide  matrix 
composite  N24A  Ml  SC/SC,  figure  5.  The  composite  used  in  this  study  has  eight  plies  for  a 
minimum  thickness  of  .08  in. 


Figure5:  Photomicrograph  of  N24A  Ml  SC/SiC  composite.  [Mital  2009] 


The  fiber  tows  are  bundles  of  Syl ramie1  -iBN  fiber  which  are  then  woven  into  a  five-harness 
fabric  perform,  figure  6a,  with  a  I  ay  up  sequence  of  [0°/90°]4s  •  There  are  20  ends  per  inch  (epi). 
Thefiber  tow  isforged  in  an  elliptical  shape  with  atotal  thickness  of  0.005  in  and  an  aspect  ratio 
of  9:1. 

Fiber 


—  BN  coating 


SiC  +  porosity 

(a)  (b) 

Figure  6:  (a)  five-harness  woven  fabric  perform,  (b)  composition  of  afiber  tow  [Mital 
2009] 


/ 

/ 
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The  fiber  tow  is  itself  a  composite,  figure  6b,  comprised  of  bundles  of  Sylramic1  -iBN  filaments 
coated  with  a  silicon-doped  boron  nitride  (BN)  interface  The  preform  isthen  impregnated  with  a 
SC  matrix  via  a  chemical  vapor  infiltration  (CVI)  process  to  complete  the  rigid  perform  mat. 
These  tows  can  be  seen  in  Figure  5  as  the  dark  grey  material.  The  preform,  Figure  6a,  is  then 
filled  with  a  SC  matrix  utilizing  a  slurry-casting  melt  infiltration  process,  which  is  visible  in 
figure  5  as  the  light  grey  material.  The  analysis  tool  may  be  used  to  model  other  material 
systems 

3.3.0  Components  of  methodology/tool 

This  methodology/tool  uses  preexisting  specialty  software  packages  which  have  been  extensively 
used  by  NASA  and  private  industry.  Application  of  this  methodology  has  been  demonstrated  on 
the  CMC  of  high  interest.  First  the  tow  properties  are  determined  from  the  properties  of  its 
constituents  using  MAC/GMC.  Then  the  laminate  properties  are  generated  using  pcGINA  or 
MAC/GMC  depending  on  the  users’  preferenca  Finally,  the  material  properties  generated  from 
pcGINA  or  MAC/GMC  input  into  the  Peri  dynamics  model  and  the  impact  simulation  is 
conducted;  which  generates  a  damage  map  and  an  estimatefor  residual  strength. 

3.3.1  Peri  dynamic  theory 

Despite  the  da/el  opment  of  many  important  concepts  to  predict  material  behavior  and  failure,  the 
prediction  of  failure  modes  and  residual  strength  of  composite  materials  is  a  challenge  within  the 
framework  of  the  finite  element  method  (FEM).  Furthermore,  the  previous  methods  cannot 
address  the  nudeation  of  damage  in  a  continuous  material.  The  field  of  fracture  mechanics  is 
primarily  concerned  with  the  evolution  of  pre-existing  defects  within  a  body,  rather  than  the 
nudeation  of  new  defects  Even  when  addressing  crack  growth,  the  FEM  suffers  from  the 
inherent  limitation  that  it  requires  remeshing  after  each  incremental  crack  growth.  In  addition  to 
the  need  to  remesh,  existing  methods  for  fracture  modding  also  suffer  from  the  requirement  of 
an  external  crack  growth  criterion.  This  criterion  prescribes  how  damage  evolves  a  priori  based 
on  local  conditions,  and  guides  the  analysis  as  to  when  and  how  damage  initiates  and  propagates 
Considering  the  difficulty  in  obtaining  and  generalizing  experimental  fracture  data,  providing 
such  a  criterion  for  damage  growth,  espedally  in  composite  structures  dearly  presents  a  major 
obstade  to  fracture  modding  using  conventional  methods  This  prevents  such  methods  from 
bang  applicable  to  problems  in  which  multiple  damage  growth  occurs  and  interacts  in  complex 
patterns 

The  governing  equations  of  the  FEM  are  based  on  the  partial  differential  equations  (PDEs)  of 
dassical  continuum  mechanics  and  that  the  spatial  derivatives  required  by  the  PDEs  do  not,  by 
ddinition,  exist  at  crack  tips  or  along  crack  surfaces  Therdore,  the  basic  mathematical  strudure 
of  the  formulation  breaks  down  whenever  a  crack  appears  in  a  body.  Various  special  techniques 
have  been  devdoped  in  fracture  mechanics  to  deal  with  this  limitation.  Generally,  these 
techniques  involve  redefining  a  body  in  such  a  way  as  to  exdude  the  crack,  then  applying 
condi  ti  ons  at  the  crack  surf  aces  as  boundary  condi  ti  ons 

In  order  to  overcome  these  problems  [Slling  2000]  introduced  a  nonlocal  theory  that  does  not 
require  spatid  derivatives — the  peridynamic  (PD)  theory.  The  man  difference  between  the  PD 
theory  and  dassical  continuum  mechanics  is  that  the  former  is  formulated  using  integral 
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equations  as  opposed  to  derivatives  of  the  displacement  components.  Thisfeeture  allows  damage 
initidion  and  propagation  at  multiple  s'tes  with  arbitrary  paths  inside  the  material,  without 
resorting  to  special  crack  growth  criteria  In  the  PD  theory,  internal  forces  are  expressed  through 
nonl  ocal  i  nteracti  ons  between  the  materi  al  poi  nts  withi  n  a  conti  nuous  body,  and  damage  i  s  a  part 
of  the  constitutive  model.  Interfaces  between  dissimilar  material  shave  their  own  properties  and 
damage  can  propagate  when  and  where  it  is  energetically  favorablefor  it  to  do  so. 

With  the  PD  theory,  damage  in  the  material  is  simulated  in  a  much  more  realistic  manner 
compared  to  the  classical  continuum-based  methods.  The  broken  interactions  may  align 
themselves  along  surfaces  that  form  cracks  and  the  deformation  is  discontinuous  across  such  as  a 
crack,  yet  the  integral  equations  continue  to  remain  valid.  The  PD  theory  has  been  utilized 
successfully  for  damage  prediction  in  many  problems. 

In  a  peri  dynamic  simuldion,  the  distances  between  points  in  the  material  are  monitored  and 
stretches  between  these  points  are  determined.  When  the  stretch  between  a  pair  of  points 
exceeds  the  specified  failure  stretch,  the  peridynamic  bond  between  these  material  points  is 
broken.  Hence,  these  material  points  will  no  longer  interact  and  failure  takes  place  Material 
failure  is  therefore  intrinsic  to  the  theory,  and  there  is  no  need  for  any  complex  treatment  to 
capture  crack  initiation  and  propagation.  Thus,  damage  is  simulated  in  a  much  more  realistic 
manner  compared  to  class  cal  continuum-based  methods 

In  the  case  of  an  isotropic  material  where  there  is  no  directional  dependence  of  the  interactions 
between  the  material  points  However,  if  anisotropic  materials  such  as  a  fiber-re'nforced 
composite  structure  are  considered,  the  directional  dependency  must  be  included  in  the 
peridynamic  analysis  Therefore,  four  different  peridynamic  material  parameters  are  introduced, 
as  shown  in  figure  7,  to  model  a  fiber-reinforced  composite  laminate  A  material  point  of  interest 
can  only  interact  with  a  material  point  located  either  in  the  same  ply  or  in  the  adjacent  plies 


Figur  e  7:  Peri  dynami  c  bonds  for  a  f  i  ber-rei  nf  orced  composite  I  ami  nate 
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3.3.2  MAC/GMC 

MAC/GMC  is  a  software  package  for  the  analysis  of  composite  material  and  laminates 
developed  at  the  NASA  Glenn  Research  Center.  It  is  based  on  the  generalized  method  of  cells 
(GMC)  micromechanics  theory  [Paley  and  Aboudi,  1992],  It  calculates  the  local  stress  and  strain 
fields  generated  in  the  composite  material  given  spatially  varying  constituents  With  this 
software  package  one  may  assess  the  inelastic  materi  al  behavi  or  and  vari  ous  type  of  damage  to  a 
laminate  given  inhomogeneous  local  models  of  the  materials  constituents  The  set  of  tools 
encoded  in  MAC/GMC  has  se/eral  advantages  over  the  commonly  U9ed  finite  element 
mi  cromechanics  approaches  MAC/GMC  provides  sa/eral  constitutive  models  for  heterogeneous 
materials  This  software  package  allows  one  to  represent  the  composite  as  continuum  without 
any  boundaries  Thus  material  models  generated  in  MAC/GMC  can  be  seamlessly  embedded 
within  higher  scale  structural  analyses  such  as  finite  element.  The  MAC/GMC  model  is 
designed  to  represent  the  composite  behavi or  at  amaterial  point. 


z 

i 


Layer 

3 

*3 
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GMC  RUC 
X 


0  Integration  Points  (+) 

(b) 


Figure8:  a)  Defining  a  Repeeting  Unit  Cell  (RUC) ;  b)  Geometry  employed  in  MAC/GMC  for 
laminate  theory  analysis  Reprinted  from  [Bednarcyk  2002] 

With  the  repeating  unit  cell  (RUC)  method,  one  identifies  a  pattern  which  repeats  throughout  the 
material.  In  the  doubly  periodic  (RUC)  method,  as  depicted  above,  the  unit  cell  repeats  in  a 
plane  whereas  with  the  triple  periodic  RUC  method  the  unit  cell  repeats  in  three  dimensions 
There  may  be  more  than  one  distinct  repeeting  unit  cel  I  in  the  material.  The  (RUC)  is  not  meant 
to  model  alaminate,  but  instead  it  considers  an  idealized  compositein  an  infinite  plane  or  space 
Thus  when  using  the  RUC  method  in  the  MAC/GMC  code,  one  is  not  concerned  with 
thicknesses  heights  and  widths 


Inside  the  RUC  of  figure  8a),  one  can  see  tows  in  both  the  fill  and  warp  directions  as  well  as 
space  where  the  matrix  material  is  otherwise  present.  Thus  the  RUC  is  further  divided  into  sub 
cells  such  that  there  is  to  be  only  one  type  of  material  per  sub  cell.  The  micromechanics 
equati  ons  are  then  appl  i  ed  and  sol  ved  for  each  sub  cel  I  wi  thi  n  the  RUC. 
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In  addition  one  may  use  MAC/GMC  to  generate  the  material  response  of  a  ply  or  laminate  via 
classical  lamination  theory,  figure  8b).  The  input  for  the  laminate  solver  is  very  similar  to  the 
input  for  the  RUC  method.  The  main  difference  being,  that  the  user  must  supply  information 
regarding  the  orientations  of  the  fibers,  thicknesses  and  number  of  plies  which  comprise  the 
laminate  pi  ate 

Previous  work  has  shown  that  the  traditional  one  step  three-di mens' onal  homogenization 
procedure  for  modeling  woven  composites  is  inaccurate  due  to  a  lack  of  shear  coupling  in  the 
model.  This  initial  shortcoming  has  been  resolved  by  the  introduction  of  the  two  step 
homogenization  procedure  Wherein,  each  RUC  identified,  is  first  homogenized  through  its 
thickness  before  being  reassembled  with  the  other  RUCs  anisotropic  D  matrices  and  then  all  are 
homogeni  zed  across  the  pi  ane 


*3 

(a)  (b) 

Figure  9:  a)  apian  weave  composite  b)  the  exploded  view  of  that  plain  weave  compos  te 
reprinted  from  [Bednarcyk  2000] 


Figure9a),  is  a  MAC/GMC  model  of  apian  weave  CMC,  figure9b)  is  an  exploded  view  of  the 
same  model.  Imagine  dividing  the  material  into  serial  cross-sections  and  separating  it  as  in 
figure9  b).  Closer  inspection  of  figure  9b)  allows  one  to  identify  six  basic  groups  which  can  be 
seen  in  figure  10  a).  A  first-step  model  is  developed  for  each  group,  and  once  MAC/GMC  has 
derived  the  full  anisotropic  D  matrix  for  each  of  the  six  groups  then  one  uses  that  information  as 
i  nput  i  nto  the  second-step  model ,  whi  ch  can  be  seen  i  n  f  i  gure  1 0b) 


Group  4 


Group  5 


Group  6 


Group  1 


Group  2  Group  3 

fa) 


(b) 


Fi  gure  10:  aO  six  distinct  materi  al  groups  f  a  the  compos  te  b)  the  reassembl  y  of  the  groups  i  n  1 0a) 
repri  nted  f  rom  [  Bednarcyk  2000] 
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3.3.3  pcGINA 

pcGINA  stands  for  “  PC  based  Graphical  Integrated  Numerical  Analysis  program”.  It  is  used  for 
designing  and  predicting  the  properties  of  textile  reinforced  composites.  It  predicts  the  static 
elastic  and  thermal  properties  of  textile  composites.  It  was  selected  for  use  in  this  work  due  to  the 
fact  that  it  is  widely  known  in  industry.  As  its  use  is  a  good  fall -back  position  and  double  check 
of  the  GM  C/M  AC  models. 

However,  it  was  developed  for  Polymer  Matrix  composites  (PMC)  not  CMCs.  It  is  not  a  micro¬ 
mechanics  code  like  GM  C/M  AC  rather  it  is  a  fiber  weave,  or  textilecode.  Thus,  there  areafe/v 
times  when  the  user  has  to  “trick”  the  user  input  screens;  as  a  PMC  does  not  have  the  extra  step 
of  a  CVI  infiltration  of  the  yarn  to  create  a  tow.  Some  DOD-AF  [Gowayed  et  al  2007]  work  has 
adapted  pcGINA  to  work  with  CMC  materials,  however  the  authors  are  not  aware  of  the  full 
extent.  None-the-less,  it  clearly  belongs  in  a  suite  of  codes  along-side  GMC/MAC,  for 
verification  and  validation  purposes. 

3.4.0  Models 

First  to  be  described  is  the  MAC/GMC  model  used  to  determine  the  tow  properties.  Then  the 
models  used  to  determine  the  laminate  material  properties  presented.  Lastly,  the  impact  models 
are  described. 

3.4.1  Tow  model 

Unlike  a  polymer  matrix  composite  (PMC),  one  must  derive  the  material  properties  of  the  tow 
before  one  may  derive  the  material  properties  of  a  ply  or  laminate  panel.  Snce,  the  tow  is 
comprised  of  filaments,  possibly  braided  into  a  yarn,  which  are  coated  by  a  boron  nitride  (BN) 
interface  material  before  bang  infiltrated  with  a  silicon  carbide  (SC)  matrix  material.  The  tow 
itself  should  bethought  of  as  a  unidirectional  lamina  So  a  dear  distinction  must  be  maintained 
between  the  ratio  of  filamentsto  chemical  vapor  infiltrated  (CVI)  matrix  material  within  the  tow, 
and  the  ratio  of  tows  to  melt  infiltrated  matrix  material  in  the  ply  or  laminata 

Filament  volume 

fraction  within  _  (number  of  filaments  )  *  (area  of  a  filament)  ^ 

the  tow  ~ 

area  of  tow 

Even  though  one  should  think  of  the  CMC  tow  as  a  unidirectional  lamina,  in  reality  it  may  not 
be,  since  the  filaments  may  be  braided  into  a  yarn.  If  this  happens  to  be  the  physical  situation 
MAC/GMC  has  models  in  its  standard  library  which  can  represent  a  braided  yarn  structure. 
Figurell  (b)  is  the  MAC/GMC  model  used  to  describe  the  idealization  of  the  tow  as  depicted  in 
figurell(a).  One  uses  equation  1  to  determine  the  filanent  volumefraction  within  the  tow. 

The  input  deck  for  the  tow  model,  figure  11b: 

Number  of  materials:  3 

Constitutive  models:  Elastic 

Materials:  iBN-Syl  ramie 

CVI-Sc 
CVI -BN 
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Analysis  type: 

Analysis  model: 

RUC  architecture: 

Fi  ber  vol  ume  f  racti  on : 
Ratio  of  interface 
thi  ckness  to  f  i  ber  radi  us: 


(*RUC)  Repeating  Unit  Cell  Analysis 
Doubly  periodic  GMC 

7x7  Circular  Fiber  Approx.,  Rectangular  Pack 
0.5 

0.124 


Figure'll:  a)  idealized  model  of  atow,  b)  MAG/GMC  model  of  the  idealized  tow 


Table  1:  Input  valuestothetow  model,  constituent  properties 


ea 

(GPS) 

Et 

(GPa) 

Va 

vr 

Ga 

(GPa) 

a 

10’6K 

Material  properties  at 

room 

temperature 

Syl  ramie 
fiber 

380 

380 

0.17 

0.17 

162.4 

4.6 

CIV-SC 

matrix 

380 

380 

0.17 

0.17 

162.4 

4.6 

BN  coating 

21 

21 

0.22 

0.22 

8.607 

5.2 

Material 

properties  at  hig 

temperature  (1 204  C) 

Syl  ramie 
fiber 

365 

365 

0.17 

0.17 

156 

8.0 

CVI-SC 

matrix 

358 

358 

0.17 

0.17 

153 

9.0 

BN  coating 

14 

14 

0.22 

0.22 

5.738 

10.0 

Table  1  lists  the  input  values  for  thi  stow  model,  figure  11b).  EA  is  Young’s  modulus  in  the  axial 
direction  while  Er  is  Young’s  modulus  in  the  transverse  direction.  Similarly,  vA  Poisson’s  ratio  in 
the  axial  direction  and  vt  is  Poisson’s  ratio  along  the  transverse  direction.  GA  is  the  shear 
modulus  in  the  xy  plane,  a  isthe  coefficient  of  thermal  expansion. 
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3.4.2  pcGINA  laminate  model 

If  the  laminate  is  balanced,  then  the  material  properties  of  the  ply  are  the  same  as  the  material 
properties  of  the  laminate,  with  the  exception  of  the  effect  of  thickness 


Figure  12:  pcGINA  model  of  a  5-harness  satin  weave  composite 


In  figure  12,  for  both  the  red  and  blue  tows,  it  is  assumed  that  the  major  axis  is  45  milli -inches 
(mil)  and  that  the  minor  axis  is  5  (mil).  The  tow  initially  starts  with  a  cylindrical  geometry,  but 
ends  up  wi  th  an  el  I  i  psoi  dal  shape  due  to  fabri  cati  on.  The  angl  e  between  the  tows  i  s  90°.  Thi  s 
leads  to  a  total  fiber  volume  fraction  of  0.70911,  of  which  0.35455  is  the  blue  tows,  and  the  other 
0.35455  is  the  red  tows  In  this  context,  fiber  does  not  refer  to  the  filaments  which  comprise  the 
tows,  instead  the  term  fiber  refers  to  the  tow.  Table  2  lists  the  input  values  used  in  the  pcGINA 
ply  model  depicted  in  Figure  12. 


Table 2:  input  material  proper lies 


Ea 

(GPa) 

Et 

(GPa) 

VA 

VT 

Ga 

(GPa) 

OA 

10‘6K 

OlT 

10'6K 

Material  properties  at 

room  temperature 

Tow  fiber 

320 

149 

0.176 

0.15 

59.6 

4.6 

4.7 

MI-SC 

matrix 

310 

310 

0.17 

0.17 

132.48 

4.7 

4.7 

aterial  properties  at  hi g 

temperature  (1204  C) 

Tow  fiber 

306 

120 

0.175 

0.137 

47.6 

8.5 

8.8 

MI-SC 

matrix 

276 

276 

0.17 

0.17 

117.95 

9 

9 

The  version  of  pcGINA  used  in  this  study,  v.11.98,  is  designed  to  model  PMCs  as  opposed  to 
CMCs  Therefore,  certa'n  manipulations  are  required:  ^ter  inputting  values  for  the  filaments, 
pcGINA  provides  a  screen  to  input  yarn  specific  parameters  The  values  inputted  on  the  yarn 
screen  are  required  by  pcGINA,  however  their  engineering  value  has  been  made  redundant  by 
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the  tow  model.  It  is  important,  that  the  values  on  this  screen  be  inputted  in  such  away  that  the 
filament  packing  factor,  generated  by  pcGINA,  is  calculated  as  1.  The  program  will  issue  a 
warning,  but  this  is  to  be  ignored.  Recall  that  the  tow’s  properties  were  calculated  prior  to  the 
use  of  pcGINA. 

3.4.3  MAC/GMC  laminate  model 

MAC/GMC  was  used  to  analyze  the  woven  reinforcement  architecture  of  theN24A  CMC.  Both 
the  one-step  and  two-step  homogenization  procedures  from  [Mital  et  al  2009]  have  been  re¬ 
produced.  The  two  step  procedure  was  implemented  since  the  single-step  procedure  has  been 
shown  to  derive  less  than  accurate  results  for  woven  composites  Conversely,  for  unidirectional 
composites  the  one-step  homogenization  procedure  provides  accurate  results  Table  2,  above, 
lists  the  material  properties  used  in  these  models  MAC/GMC  internal  constitutive  model  sal  low 
specification  of  an  arbitrary  direction  of  tra/erse  isotropy.  Material  5  is  the  melt  infiltrated  SC 
matrix.  Materials  1-4  and  6-7  are  all  the  fiber  tow  and  thus  have  the  same  material  properties 
The  di  fference  between  them  i  s  a  vector  whi  ch  defines  thei  r  di  recti  on  of  traverse  i  sotropy . 

The  i  nput  deck  for  the  woven  CM  C  model .  f  i  gure  1 4: 

Number  of  materials  7  (NMATS=7) 

Constitutive  models  Arbitrary  transversely  isotropic  (CMOD=9) 

Elastic  (CMOD=6) 

Materials  User-defined  (M  ATI  D=U) 

Material  property  source:  Read  from  input  file  (MATDB=1) 

Material  properties  See  table  2  above 

Direction  of  trans  Isotropy: 

Material  #1:  (0,  0, 1) 

Material  #2:  (0, 1,0) 

Material  #3:  (2.5, 14,  0) 

Material  #4:  (-2.5, 14,  0) 

Material  #6:  (2.5,  0, 14) 

Material  #7:  (-2.5,  0, 14) 

(*  RUC)  Repeating  Unit  Cell  Analysis 
Triply  periodic  GMC  (MOD=3) 

User-defined  (ARCHID=99) 

4  (NA=4) 

10  (NB=4) 

10  (NG=10) 

2.5,  2.5,  2.5,  2.5  (D=2.5,  2.5,  2.5,  2.5) 

36.0,1 4.0,36.0, 1 4.0,36.0, 1 4.0,36.0, 1 4.0,36.0, 1 4.0 
36.0, 1 4.0,36.0, 1 4.0,36.0, 1 4.0,36.0, 1 4.0,36.0, 1 4.0 


Analysis  type: 

Analysis  model: 

RUC  architecture: 

No.  sub-cellsin  xl-dir.: 
No.  sub-cellsin  x2-dir.: 
No.  sub-cellsin  x3-dir.: 
Sub-cell  depths: 
Sub-cell  heights 
Sub-cell  lengths 
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a)  b) 

Figure  13:  a)  material  coordinate  system  b)  legend  for  fiber  orientation 

Fi  gures  1 3a)  and  1 3b)  are  the  materi  al  coordi  nate  system  and  the  I  egend  f  or  the  f  i  ber  ori  entati  on 
of  figure  14.  The  vector  coordinates  listed  in  the  legend  13b)  are  defined  by  the  coordinate 
system  depicted  in  figure  13a).  The  fibers  where  not  draw  on  the  black  and  purplecells;  as  these 
lines  would  have  to  be  drawn  through  the  thickness  of  the  page.  The  green  cells  represent  the 
melt  infiltrated  matrix  material.  The  red,  orange  and  yellow  cells  represent  the  tows  in  the  x2 
direction,  whereas  the  blue,  black  and  purple  cel  Is  are  the  tows  in  thex3  direction.  Inspection  of 
gamma  one  in  figure  14,  starting  at  the  leftmost  column  in  thex2  direction,  shows  the  red  tow  at 
the  bottom  two  rows  As  one  moves  i  n  the  x2  di  recti  on,  the  red  tow  stays  i  n  the  bottom  two  rows, 
until  the  sixth  column  is  reached,  at  which  point  the  red  fiber  tow  begins  to  rise,  depicted  as  the 
orange  cells,  and  then  in  the  seventh  column  the  red  fiber  tow  is  on  the  topmost  two  cells  in  the 
xl  direction.  In  the  eighth  column  along  thex2  direction,  the  red  fiber  tow  begins  to  drop  down 
in  thexl  direction,  which  is  depicted  as  the  yellow  cells  Finally  the  red  tow  is  once  again  on  the 
bottom  most  rows  for  columns  nine  and  ten.  The  same  behavior  can  be  seen  for  the  blue,  black 
and  purple  tow.  However,  instead  of  staying  on  one  gamma  cross-sectional  slice,  the  weaving 
behavior  of  the  blue  tow  will  be  seen  by  traversing  from  gamma  one  to  gamma  five  Consider 
gamma  one,  column  three,  the  top  two  rows  in  this  location  is  the  blue  tow.  Now,  in  gamma 
two,  the  bl  ue  tow  i  s  depi  cted  as  purpl  e  cel  I  s  as  i  t  drops  down  i  n  the  x  1  di  recti  on.  Then  i  n  gamma 
three,  the  blue  tow  occupies  the  two  bottommost  cells  in  the  xl  direction.  In  gamma  four,  the 
blue  tow  rises  through  thexl  direction  and  this  is  depicted  by  the  black  cells  Finally,  in  gamma 
five,  the  blue  tow  has  returned  to  the  topmost  cells  in  thexl  direction. 

Thus  one  may  build  a  MAC/GMC  model  represent  any  custom  material  architecture  The  input 
deck  has  commands  to  define  each  material  of  interest.  One  then  has  a  set  of  commands  to  define 
how  the  unit  cells  are  connected  and  which  material  they  represent.  Finally,  one  has  control  over 
the  dimensions  of  each  unit  cell.  Recall,  that  the  MAC/GMC  model  represents  a  material  point 
which  is  callable  by  a  structural  model.  Embedded  in  MAC/CMC  is  a  library  of  material 
constitutive  behavior  models  This  eliminates,  in  many  cases,  the  need  to  write  custom 
subroutines 
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Rgure12:  MAC/GMC  model  of  a5harnesssatinwea/eCMC 
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3.4.4  PD  model  for  a  FOD  event: 

A  9' C/SC  ceramic  matrix  composite  (CMC)  is  impacted  either  by  a  rigid  or  a  deformable  (steel) 
indenter  as  shown  in  figure  13.  Two  different  boundary  conditions  are  considered,  i.e.  full 
support  and  partial  support. 

The  geometrical  properties  of  the  SC/SC  CMC  is  specified  as  L  =  45mm,  W  =  8mm , 
H  =  2.25  mm.  The  spherical  impact  or  has  a  diameter  of  D  =  1.5mm.  For  the  partial  support 
case,  the  distance  between  two  supports  is  specified  as  B  =  20mm. 

The  SC/SC  CMC  material  is  assumed  to  be  isotropic  with  an  elastic  modulus  of  E  =  220GRa 
and  density  of  /?  =  2400kg/m3.  The  elastic  modulus  and  density  of  the  steel  impactor  are 
E  =  200 GRa and p  =  7860 kg/m3,  respectively.  The  steel  impactor  can  exhibit  plastic 
deformation  when  it  experiences  stress  state  near  the  ultimate  strength  of  the  material  which  is 
specified  as  2.0  GRa  The  9C/SC  CMC  and  steel  impactor  have  fracture  toughness  values  of 
K,  =  25. M PaVm  and  K,  =  90 M RaVm  ,  respectively. 


(b) 


(a) 


Figure  13:  S  C/9  C  CMC  is 'impacted  by  a  rigid  indenter  (a)  Full  support  (b)  Partial  support 

Two  different  impactor  velocity  conditions,  v0  =  200 nVs  and  v0  =  400 nVs  are  considered  for 
both  full  and  partial  support  cases. 

Peri  dynamic  model  is  generated  by  using  51976  material  points  with  grid  size  and  horizon  values 
of  Ax  =  0.25mm  and  5  =  0.75375mm,  respectively.  The  bond  constants  for  the  CMC  and  steel 

impactor  are  c  =  2.6x1012  N/m6  and  c  =  2.37x1012  N/m6,  respectively.  The  critical  stretch 
values  of  CMC  and  steel  impactor  are  specified  as  ^  =  0.00276  and  =0.0157,  respectively, 
whi  ch  can  be  computed  by  usi  ng  the  rel  ati  ons  deri  ved  by  S 1 1  i  ng  and  A  skari  [8] , 

Sq=  4K2'  -  for  2-D  CMC  plate  and  ^  =  - '^i-  for  3-D  steel  impactor 
y9E  £  y3E  £ 

9able  time  step  size  used  in  the  dynamic  analysis  is  chosen  as  At  =  1.3xl0*s.  The 
displacement  boundary  conditions  are  applied  to  one  layer  of  material  points  at  the  bottom  edge 
of  the  model  by  constraining  all  displacement  components. 
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4.0  Results 


4.1.0  MAC/GMC  tow  results 

Table  3  lists  a  comparison  between  the  results  of  this  work  and  the  results  published  in  [Mital  et 
al  2009].  There  are  slight  discrepancies  between  the  two  sets  of  results.  However,  most  of  these 
differences  are  less  than  10  percent.  Ea  is  Young’s  modulus  in  the  axial  direction  while  Et  is 
Young’s  modulus  in  the  transverse  direction.  Similarly,  va  Poisson’s  ratio  in  the  axial  direction 
and  vt  is  Poisson’s  ratio  along  the  transverse  direction.  Ga  is  the  shear  modulus  in  the  xy  plane, 
a  is  the  coefficient  of  thermal  expansion. 


Table 3:  Comparison  of  theMAC/GMC  tow  model  resultswith  thepublished  results 


Ea 

(GPa) 

Et 

(GPa) 

vA 

VT 

Ga 

(GPa) 

ClA 

10'6K 

ClT 

10"6K 

Literature  resull 

sat  room  temperature 

320 

149 

0.176 

0.15 

59.6 

4.6 

4.7 

This  study’s  results  at  room  temperature 

320.8 

150.4 

0.1756 

0.1507 

60.12 

4.662 

5.328 

Literature  results  at  high  tern 

perature(1204C) 

306 

120 

0.175 

0.137 

47.6 

8.5 

8.8 

This  study’s  results  at  high  temperature  (1204  C) 

304.7 

120.5 

0.1752 

0.1378 

47.94 

8.409 

8.744 

4.2.1  pcGINA  laminate  results 

Tables  4  and  5  ha/e  a  comparison  of  resultsfrom  thisstudy  and  the  results  published  in  [Mital  et 
al  2009].  The  difference,  between  the  material  properties,  generated  by  the  two  studies  are  for 
most  case  within  rounding  error.  The  worst  discrepancy  is  Young’s  modulus  in  the  z  direction 
for  the  room  temperature  case.  The  two  differ  by  approximately  4.6%. 


Table  4:  Comparison  of  room  temperature  results wi 


h  the  literature 


Ex 

(GPa) 

Ey 

(GPa) 

Ez 

(GPa) 

Vxy 

Ga 

(GPa) 

Thisstudy 

247.604  ' 

247.604 

182.146 

0.1138 

103.635 

Lit.  model 

248 

248 

174 

0.11 

102 

Lit.  data 

252 

252 

82 

0.13 

n/a 

Table 5:  Comparison  of  high  tern 

perature(1204C)  resultswith  theliterature 

Ex 

(GPa) 

Ey 

(GPa) 

Ez 

(GPa) 

Vxy 

Ga 

(GPa) 

a 

10  6  K 

kx 

kz 

Thisstudy 

221.746  J 

221.746 

150.567 

0.116  4 

87.82 

7.811 

17  43  ! 

12.66 

Lit.  model 

221 

221 

151 

0.11 

86 

7.8 

18 

12 

Lit.  data 

230 

230 

n/a 

n/a 

n/a 

6 

20 

18 

20 


4.2.2  MAC/GMC  laminate  results 

Tabl  es  6  and  7  compare  the  resul  ts  found  i  n  thi  s  study  wi  th  the  resul  ts  publ  i  shed  i  n  I  i  terat  ure 
[Mital  et  al  2009]. 


Table6:  Sng 

le-step  homogenization  procedure 

Ex 

(GPa) 

Ey 

(GPa) 

Ez 

(GPa) 

Vxy 

Ga 

(GPa) 

Material  properties  at  room  temperature 

This  study 

232.8 

232.9 

182.8 

0.1243 

70.64 

Lit.  model 

233 

233 

183 

0.125 

70.7 

Lit.  data 

252 

252 

82 

0.13 

n/a 

M  ateri  al  properti  es  at  hi  gh  temperature  ( 1 204  C) 

This  study 

202.6 

202.7 

151.4 

0.113 

57.33  1 

Lit.  model 

203 

203 

151 

0.11 

57.3 

Lit.  data 

230 

230 

n/a 

n/a 

n/a 

T able 7:  T wo-stq)  homogenization j^ocedu re 


Ex 

(GPa) 

Ey 

(GPa) 

Ez 

(GPa) 

Vxy 

Ga 

(GPa) 

r 

Material  properties  at  room  temperature 

This  study 

252.94 

249.01 

183.11 

0.128 

74.11  ! 

Lit.  model 

252 

252 

183 

0.127 

74.5 

Lit.  data 

252 

252 

82 

0.13 

n/a 

M  ateri  al  properti  es  at  hi  gh  temperature  ( 1 204  C) 

Thisstudy 

227.78 

223.43 

151.60 

0.117 

60.89 

Lit.  model 

225 

225 

151 

0.12 

60.9 

Lit.  data 

230 

230 

n/a 

n/a 

n/a 
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4.3.1  PD  resultsfor  a  rigid  impactor 

The  damage  predictions  obtained  for  two  different  boundary  conditions  for  an  impactor  velocity 
of  v0  =  200rrVs  and  v0  =  400 rrVs are  depicted  in  Figures  14  and  15.  In  figure  15(a),  a  substantial 

damage  occurs  around  the  impacted  region  for  the  full  support  case.  For  the  partial  support  case, 
a  conical  damage  shape  is  observed  and  the  damage  occurs  not  only  at  the  impacted  region  but 
also  at  the  bottom  of  the  specimen.  For  ahigher  velocity  case,  i.av0  =400rrVs,  similar  damage 
shapes  are  observed  wi  th  a  bi  gger  damage  area 


(a)  (b) 

Figure  14:  Damage  plot  for  SC/SC  CMC  impacted  by  a  rigid  indent  er  with  a  velocity  of 
v0  =  200nVs  (a)  Full  support,  (b)  partial  support 


(a)  (b) 

Figure  15:  Damage  plot  for  SC/SC  CMC  impacted  by  a  rigid  indenter  with  a  velocity  of 
v0  =400rrVs  (a)  Full  support,  (b)  partial  support 

Peri  dynamic  damage  prediction  can  be  qualitatively  compared  against  experimental  results 
obtained  by  [Choi  2008],  as  shown  in  Figure  3.  The  experimental  observations  and  PD 
simulation  results  appear  to  be  in  dose  agreement  as  shown  in  Figures14and  15. 
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4.3.2  PD  resultsfor  a  deformable  (steel)  impact  or 

The  damage  plots  obtained  for  two  different  boundary  conditions  for  an  impactor  velocity  of 
v0  =  200 nVs  and  vo  =  4O0nVsare  depicted  in  figures  16  and  17.  As  shown  in  figure  16a, 

substantial  damage  occurs  around  the  impacted  region  for  the  full  support  case.  For  the  partial 
support  case,  in  addition  to  the  damage  occurrence  around  the  impacted  region,  a  coni  cal  damage 
occurs  at  the  bottom  of  the  sped  men  as  shown  i  n  f  i  gure  1 6b.  For  both  cases,  the  steel  i  mpactor 
is  significantly  damaged  and  experiences  plastic  deformation.  For  a  higher  velodty  case,  i.a 
vo=4O0nVs,  shown  in  figure  17,  similar  damage  pdterns  are  observed  with  larger  damage 


Figure  16:  Damage  plot  for  SC/SC  CMC  impacted  by  a  steel  impactor  with  a  velodty  of 
v0  =  200rrVs  (a)  Full  support,  (b)  partial  support 


Figure  17.  Damage  plot  for  SC/SC  CMC  impacted  by  a  steel  impactor  with  a  velodty  of 
v0  =400rrVs  (a)  Full  support,  (b)  partial  support 
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5.0  Conclusions 

There  is  substantial  interest  in  developing  CMC  material  systems  for  U9e  in  turbo-machinery 
appl i cati ons  There  are  def i ni te  i ncrea9es  to  sped f i c  power  and/or  reducti ons  in  fuel  consumpti on 
when  CM  Cs  are  i  incorporated  i  nto  aeroengi  ne  desi  gn. 

The  objective  of  phase  I  was  to  develop  and  demonstrate  an  initial  methodology  for  analyzing 
the  effects  of  f  orei  gn  obj  ect  damage  (FOD)  on  aeroengi  ne  components  uti  I  i  zi  ng  CM  C  materi  al  s 

To  engineer  FOD  resistant  materials,  one  needs  models  which  are  able  to  simulate  a  proper 
damage  map  and  to  assess  the  residual  strength  of  the  material.  To  accomplish  this  task  the 
models  need  to  properly  represent  crack  formation  and  growth.  Proper  control  of  the  fiber/matrix 
interface  characteristic  is  the  key  to  understanding  crack  propagation  and  thus  obtaining  a  tough 
ceramic  compos  te. 

This  objective  has  been  successfully  completed  by  the  development  of  a  proof-of -concept  model 
for  one  particular  CMC  material  system.  The  effective  material  properties  of  the  CMC  material 
system  were  generated  by  using  MAC/GMC  to  homogenize  the  constituent  properties.  A 
Peri  dynamic  model  was  used  to  simulate  the  FOD  impact,  and  the  results  match  literature 
ex  peri  ments  very  wel  I . 

Feridynamic  theory  is  an  evolution  of  the  finite  element  method.  Its  main  strength  over  finite 
elements  is  that  the  damage  law  is  embedded  into  a  given  model’s  “mesh”.  Thereby,  allowing 
damage  to  occur  along  the  paths  which  are  energetically  favorable.  Whereas,  in  theFEA  method, 
one  is  required  to  constantly  re-mesh  the  model  and  have  apriority  information  about  the  crack 
path. 

The  suite  of  integrated  codes  demonstrated  in  this  effort  offer  a  cost  effective  and  reliable  tool 
for  design.  The  use  of  this  tool  will  reduce  number  of  experi mental  tests  required  early  in  the 
design  processfor  an  application  built  with  CMC  materials. 

The  presentation,  at  the  2010  CMC  session  of  Cocoa  Beach,  of  phase  I  results  generated  much 
interest  from  potential  industrial  partners,  such  as  General  Electric  and  Pratt  &  Whitney 
RocketDyna 

6.0  Recommendations 

Whi  I  e  the  phase  I  resul  ts  are  encouragi  ng,  f  urther  work  i  s  requi  red  to  f  ul  I  y  devel  op  and  opti  mi  ze 
this  approach.  The  next  obvious  step  is  to  demonstrate  the  approach  using  pertinent  data  obtained 
f  rom  vari  ous  materi  al  s  systems  together  wi  th  a  range  of  i  mpact  vari  abl  es. 

To  ensure  commercial  relevance  we  intend  to  develop  a  strategic  partnership  with  GE  or 
RocketDyne;  whi  ch  wi  1 1  guarantee  that  the  methodol  ogy  i  s  appl  i  ed  to  aero-engi  ne  appl  i  cati  ons  of 
interest  to  NAVAIR. 

The  phase  I  literature  survey  indicated:  the  need  to  explicitly  model  the  effect  that  EBC  coatings 
have  with  respect  to  FOD;  the  need  to  explicitly  model  the  effects  of  voids,  present  in  the  matrix 
materials 
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As  this  methodology  matures  there  will  beaneed  to  increase  the  coupling  of  MAC/GMCandPD 
and  prototype  the  automation  of  data  transfer  between  the  two  software  packages. 

7.0  Short  abstract  of  proposed  follow-on  research  and  development 
The  priority  objective  is  to  facilitate  the  development  of  CMC  material  systems  for  use  in  turbo¬ 
machinery  applications  In  support  of  this  goal,  the  obvious  step  is  to  analyze  an  aeroengine 
component  with  the  phase  I  methodology.  Thus  a  component  of  interest  will  be  decomposed 
into  a  set  of  streamlined  geometries  that  are  solved  piecewise  by  the  phase  I  tool.  A  detailed 
analysis  us  ng  Peri  dynamics  models  of  various  CMC  material  systems  would  be  conducted.  This 
requires  conducting  experimental  coupon  testing  to  calibrate  the  PD  impact  models  to  the 
various  material  systems  of  interest.  Then  the  optimal  material  configuration  would  be  evaluated 
in  an  appropriate  structural  model.  This  approach  ensures  that  meaningful  information  is 
generated  as  the  tool  is  refined.  If  the  loading  spectrum  is  available  then  the  fatigue  life  can  be 
estimated  us'ng  GMC/MAC. 

Many  of  the  composite  material  systems  are  novel  and  lack  extensive  databases.  Appropriate 
design  margin  can  be  established  by  probabilistic  structural  and  sens' tivity  analyses.  These 
techniques  will  also  quantify  the  effects  of  uncertainties  and  identify  influential  variables.  This 
will  ensure  that  the  material  system  is  insens  tive  to  key  design  criteria 

These  activities  will  ensure  that  the  requirements  for  a  prototype  software  process  will  be 
specified  to  control  the  automation  of  dataflow  between  software  codes. 

8.0Transtion  plan 

We  plan  to  develop  a  strategic  partnership  with  either  GE  or  RocketDyne  to  apply  the 
methodology  to  an  aeroengine  application  of  interest  to  NAVAIR.  The  phase  II  work  will 
require  a  balanced  development  effort  between  generating  a  polished  and  refined  /  re-useable 
analysis  tool  and  applying  the  current  tool  to  a  more  complex  application  of  interest  to  the 
stakeholders 

Given  an  application  problem  of  interest,  for  the  industrial  partner  and  NAVAIR,  N&R 
Engineering  and  UofA  would  work  on  both  the  tool  development  and  an  idealized  structural 
model  which  would  be  delivered  to  both  the  industrial  partner  and  NAVAIR.  The  industrial 
partner  would  be  free  to  customize  this  model  as  suits  their  needs  with  N&  R  Engineering  and 
UofA  providing  tool  support/development  as  required. 

A  marketing  study  will  be  conducted  to  evaluate  potential  interest  for  this  tool  in  other  non- 
aeroengi  ne  appl  i  cati  ons 
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200  time  step 


400  time  step 


600  time  step 


Appendix  A:  Timedependent  PD  results 


100  time  step 


300  time  step 


500  time  step 


28 


700  time  step 


800  time  step 


900  time  step 
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